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The electron-electron correlations in a many-electron [Ne = 1,2, ...,5) quantum dot confined by 
a parabolic potential is investigated in the presence of a single magnetic ion and a perpendicular 
magnetic field. We obtained the energy spectrum and calculated the addition energy which exhibits 
cusps as function of the magnetic field. The vortex properties of the many-particle wave function 
of the ground state are studied and for large magnetic fields are related to composite fermions. The 
position of the impurity influences strongly the spin pair correlation function when the external 
field is large. In small applied magnetic field, the spin exchange energy together with the Zeeman 
terms leads to a ferromagnetic-antiferromagnetic(FM-AFM) transition. When the magnetic ion is 
shifted away from the center of the quantum dot a remarkable re-entrant AFM-FM-AFM transition 
is found as function of the strength of the Coulomb interaction. Thermodynamic quantities as the 
heat capacity, the magnetization, and the susceptibility are also studied. Cusps in the energy levels 
show up as peaks in the heat capacity and the susceptibility. 
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I. INTRODUCTION. 

Magnetic-doped quantum dotsi have attracted consid- 
erable theoretical and experimental interests over the last 
two decades. Diluted magnetic II- VI and III-V semi- 
conductor (DMS) quantum dots were fabricated with 
single-electron control^i^. A rich variety of different mag- 
netic and optical properties were discovered^'^'^'^'^. 
In such structures one can explore the physical proper- 
ties coming from inter-carrier interactions and the inter- 
action of the carriers with the magnetic ion. This system 
promises to be relevant for future quantum computing 
devices, where for instance the spin of the magnetic ion 
is used as a quantum bit. More recently, electrically ac- 
tive devices were fabricated in which a single manganese 
ion is inserted into a single quantum dot^^ with control 
of the amount of charge in the quantum dot and con- 
sequently the possibility of the control of magnetism of 
single Mn-doped quantum dots. 

Investigation of the exact electronic structure of a two- 
dimensional quantum dot confined by a parabolic poten- 
tial containing several electrons and a single magnetic im- 
purity (in this paper Mit?~^) in the presence of an applied 
magnetic field is a new topic. In a recent investigationH^ 
a three-dimensional (3D) Cd(Mn)Se quantum dot con- 
taining several electrons, where only the low-energy lev- 
els of the single-electron problem were taken into account, 
was investigated in the presence of a magnetic field. Here 
we will extend this work and include all relevant energy 
levels in order to obtain a convergent solution for the 
ground state (and also excited spectrum) of the system. 

It is known that in the absence of a magnetic ion, an ex- 
ternal magnetic field is able to change the spin polarized 
state of weakly interacting electrons in a quantum dot 
in such a way that in the ground state it maximizes the 
total spin of the system: i.e. S = iVe/2. If the inter- 
particle interaction is strong, even without an applied 



magnetic field the electrons may already be polarized. 
However, with increasing magnetic field and in the case 
the inter-particle interaction is strong, the total spin of 
the system can be unusually reduced by the magnetic 
fieldl^. 

In the present study, we investigate theoretically the 
few-electron two-dimensional confined quantum dot sys- 
tem that contains a single magnetic ion in the presence of 
an external magnetic field taking into account a sufficient 
large number of single-particle orbitals such that numer- 
ical "exact" results are obtained. We explore how sensi- 
tive the whole system is to the position of the magnetic 
ion in the quantum dot and to the presence of a magnetic 
field and investigate the competition between the follow- 
ing three energies: i) the interaction of the magnetic ion 
with the electrons, ii) the interaction of the magnetic ion 
with the magnetic field, and iii) the interaction of the 
external field with the electrons. These terms affect the 
spin polarization of the N,, electrons in the quantum dot. 

Explicit studies of a iVe-correlated-electron system in- 
teracting with a single magnetic ion in nonzero mag- 
netic field are very rare in the literature. Recent 
theoreticaliSiii^'i^ii^ and experimental workiii has focused 
either on a small number of electrons using the exact di- 
agonalization approach at zero field (Refs.^W§) for a 2D 
quantum dot or at nonzero field^^ including only the low- 
est single-particle states for a 3D system or the exciton 
states relevant for optical spectroscopy of self-assembled 
magnetic-doped quantum dot a^'^d^ . 

Here, we will examine thoroughly the exact properties 
of the system containing several correlated electrons and 
a single magnetic impurity in the presence of a magnetic 
field. In our numerical "exact" diagonalization approach 
we include an arbitrary number of single-particle states to 
guarantee the accuracy of our results. We investigate the 
influence of the strength of the inter-particle interaction 
and the position of the magnetic ion on the ground state 
of the system. We predict the interesting phenomenon 
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that the magnetic ion ferromagnetically couples with the 
electrons in a region below a critical magnetic field and 
antiferromagnetically with the electrons above this crit- 
ical field. Thermodynamic properties as magnetization, 
susceptibility, and the heat capacity are investigated as 
function of magnetic field and temperature. 

This paper is organized as follows. Section II intro- 
duces the model and the numerical method. In sec- 
tion III, we present our numerical results for the many- 
particle ground state and investigate correlations through 
the appearance of vortices in the many-electron wave 
function. Sect. IV addresses the many-particle spectrum 
and in Sect. V we present results for different thermo- 
dynamic quantities. Our discussion and conclusions are 
presented in Sect. VI. 



FIG. 1: (Color online) Total Zeeman energy for a quantum 
dot with diflterent number of electrons. The magnetic ion is 
located at (0.41Zo,0), and the Coulomb interaction strength 
Ac = 0.2. The inset is the first derivative of Az with respect 
to Q,c that highlights the non-linearity of Az in certain regions 
of the magnetic field. 



II. THEORETICAL MODEL 

A quantum dot containing N,, electrons with spins 
confined by a parabolic potential and interacting with a 
single magnetic ion (Mri^+) with spin AI and a magnetic 
field is described by the following Hamiltonian: 
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The vector potential A is taken in the symmetric gauge: 
A = B/2{—y,x,0) where the magnetic field B points 
perpendicular to the plane of the interface. The confine- 
ment frequency uo is related to the confinement length 
by: ^0 = \/h/rn*ujo- 5e and gMn are the Lande g-factor 
of the host semiconductor and the magnetic ion, respec- 
tively. The dimensionless Coulomb strength is defined as: 
Ac = with a*g = 47reoeft^/m*e^ the efi'ective Bohr 

radius. The cyclotron frequency is: ujc = eB/m*. Lz, 
Sz, and Mz are the projections of the total angular mo- 
mentum of the electrons, their spins, and the magnetic 
ion in the direction of the magnetic field. The electrons 



and the magnetic impurity in the quantum dot interact 
via the contact exchange interaction with strength Jc. 

We use the set of parameters^^i^ that is applicable to 
Cd(Mn)Te which is a II(Mn)VI quantum dot with typical 
lateral size of about tens of nanometers. The dielectric 
constant e = 10.6, effective mass m* = O.lOGrno: o,*b — 
52.9 A, ge = -1.67, gMn = 2.02, Jc = 1.5 x lO^meVA^, 
and ^0 about tens of nanometers {hiUQ corresponding to 
tens of meV). For example, with Hluq — 51.32 meV gives 
?o = 26.45 A. 

We rewrite the Hamiltonian in second-quantized form: 
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where the first term is the single-particle energies Ei^^r for 
an electron in state i with spin a and the second term is 
the Coulomb interaction. The third term is the electron 
and magnetic ion Zeeman energy. The last sum is the 
electron- Mn interaction in which the first term describes 
the difference between the number of spin up and spin 
down electrons, and the last two terms describe the en- 
ergy gained by flipping the electron spin along side with 
flipping the spin of the magnetic ion. Mz, M+, and 
are the z-projection, raising and lowering operators, re- 
spectively, of the magnetic ion spin (we consider Mn ions 
which have a spin of size M=5/2). 

The single-particle states in a parabolic confinement 
potential define a complete basis of Fock-Darwin orbitals 
(/>„; (T*) and spin functions Xa 



with the Fock-Darwin orbitals: 
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In Hamiltonian i denotes a set of quantum num- 
bers {n,l} with n, I the radial and azimuthal quan- 
tum numbers, respectively. The effective length Ih — 
\J h/ra*LOH — lo/oi in the presence of a magnetic 
field is defined through the hybrid frequency ujh = 
LJo^/l + {ujc/2uJo)'^ where a = + (wc/2wo)^- The 
single-particle orbital energy is given by: 



Ei^^ = hLUH{2n + \l\ + 1) - hLUcl/2. 



(5) 



The interaction parameters between the electrons and 
the magnetic ion in the quantum dot is expressed: 
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as a product of two Fock-Darwin orbitals calculated at 
the position of the magnetic ion. 

We construct the many-particle wave function follow- 
ing the configuration interaction (CI) method. 
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where 5'^ is the k-th state of the non-interacting many 
electron wave function determined by iVg electrons with 
Ne different sets of quantum numbers {n,l,a) and the 
single scatterer with one of the six states of the Mn. 
x*={'ri, Si ) stands for the coordinates and spin of a single 
electron. In second quantization representation, the state 
^fc, which is a Slater determinant composed of single- 
electron states, can be translated into a ket vector \k) 
grouping a total of electrons into electrons with 



z-component of spin up and 
z-component of spin down: 



Np — Nt electrons with 
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FIG. 2: (Color online) Total Zeeman energy (a) calculated 
for a three-electron quantum dot with the magnetic ion at 
(0.5/o, 0) and Ac ~ 0.5. The upper inset is the first derivative 
of the main plot with respect to Qc to visualize more clearly 
the cusps. The lower inset is a zoom of the rectangular regions 
of the main plot, (b) The z-projection of the total spin of 
the magnetic impurity (solid black curve), of the electrons 
(dashed red curve), and the total angular momentum (dash- 
dotted green curve). 



Here ii| i^vt said jii j^i are the indices of the single- 
electron states for which each index is a set of two quan- 
tum numbers (radial and azimuthal quantum numbers), 
as mentioned above. 

The CI method, which is in principle exact if a suffi- 
cient number of states are included, is limited to a small 
number of electrons due to computational limitations. 
For a larger number of electrons and/or magnetic ions, 
other approaches that e.g. are based on spin-density 
functional theory (SDFT) using e.g. the local spin den- 
sity approximation (LSDA) as was used in ReL— is able 
to handle a large number of electrons and/or magnetic 
ions. The LSDA is exact only in the case of the homo- 
geneous electron gas, and in practice, works well also in 
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FIG. 3: (Color online) Zeeman energy of the magnetic ion and 
exchange interaction of a three-electron quantum dot with the 
magnetic ion located at (0.5Zo,0) and the Coulomb strength 
Ac = 0.5. The upper inset is a zoom of the small magnetic 
field region and the lower one presents the average of , Sz , 
and Lz for the ground state at very small magnetic fields. 



most inhomogeneous systems. However, in really highly 
correlated few-particle systems as discussed in this paper, 
the LSDA might fail or be at least less accurate. 



III. GROUND-STATE PROPERTIES 
A. Zeeman energy 

We first explore the magnetic field dependence of the 
total Zeeman energy: 



Az — Ec — Ejjc — ^ 
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which is the difference in energy in the presence and with- 
out a magnetic ion. It consists of three terms: the dif- 
ference of the Zeeman energy of the electrons A^'^^*'"°", 
the Zeeman energy describing the interaction of the mag- 
netic ion having spin M = 5/2 with the magnetic field, 
A^", and the exchange interaction of the ion with the 
electrons, — A|^'^. A^^"^ is just the so called "local Zee- 
man splitting term" as discussed in Ref.^^. This sum is 
basically the difference in the Zeeman energy of the elec- 
trons between the cases with and without a magnetic ion 
plus the energy contribution of the magnetic ion. 

For Nf, = 1, we find a total Zeeman energy that ap- 
pears linear in magnetic field. A similar linear behavior 
is also found in the cases with TV. > 1 but with different 
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FIG. 4: (Color online) Exchange interaction of a three- 
electron quantum dot with the magnetic ion located at 
three different positions: (0,0), (0.5io,0), and (Zo,0) for the 
Coulomb strength Ac ~ 0.5. The horizontal fine at A.^^" = 
separates the Af"'^ plane into ferromagnetic (FM) and anti- 
ferromagnetic (AFM) regions. 



slopes (see Fig. [T]). Let us suppose, in Hamiltonian ((T]), 
that the contribution from the last term (the local Zee- 
man energy or the exchange interaction term) is zero. For 
instance this is the case when a magnetic ion is located at 
the center of the quantum dot having three-electrons in 
the partially filled p shell, in which the first two-electrons 
fully fill the s shell and the remaining one is in either of 
the orbitals of the p shell. Then a perfect linear behavior 
of the total Zeeman energy is found. 

A closer look to A^ gives us a slightly different picture, 
as is provided by taking the derivative (see the inset of 
Fig. [T]). Notice that the total Zeeman term A^ has pro- 
nounced cusps and the number and positions of these 
cusps is different for different number of electrons iVg. 
There exists one at VLc = 2.6 for a two-electron quantum 
dot, one at fie = 3.4 for a three-electron quantum dot, 
and two at fie = 2.5 and 3.8 for the four-electron quan- 
tum dot, with the magnetic ion located at (0.4I^o,0). 
The three-electron system exhibits a much richer behav- 
ior when we increase the Coulomb interaction strength to 
Ac — 0.5 as seen in Fig. Ufa) where we placed the mag- 
netic ion at (0.5Zo,0). Cusps, which are highlighted in 
the two insets in Fig. [2l^a) , appear when the total angu- 
lar momentum and/or the total z-projection of the spin of 
the electrons change abruptly with magnetic field. Notice 
that the total Zeeman energy of a two-electron quantum 
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FIG. 5: (Color online) Exchange interaction of a one-electron 
(a) and two-electron (b) quantum dot with the magnetic ion 
located at three different positions and the Coulomb strength 
Ac = 0.5. The inset in the upper plot are the averages of 
the three quantities: M^, Sz, and with the ion located at 
(0.5Zo,0) at small magnetic fields. The horizontal line sep- 
arates the plane into ferromagnetic (FM) and antiferromag- 
netic (AFM) regions in both two plots. The schematic dia- 
gram in (b) explains why the exchange energy is almost zero 
in the case A'^e = 2. 



dot in the presence of the magnetic ion does not produce 
a similar behavior due to the fact that the z-projcction of 
the total spin is zero making the main contribution (from 
the Zeeman spin term of the Mn-impurity) negligible. 

The Coulomb strength and the position of the mag- 
netic impurity affects the total Zeeman energy and in- 
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FIG. 6: (Color online) Exchange interaction of a four-electron 
quantum dot with the magnetic ion located at three different 
positions and the Coulomb strength Ac ~ 0.2. The hori- 
zontal line separates the plane into ferromagnetic (FM) and 
antiferromagnetic (AFM) regions. The inset shows the total 
2-projection of the Mz, Sz, and Lz of the electron system 
with the magnetic ion located at the center of the dot. 



fluences the number and the position of the cusps. The 
first pronounced cusp appears at lower magnetic field for 
larger Coulomb interaction strength. This is a conse- 
quence of the competition between the Coulomb energy 
and the energy gap of the single-particle problem. Larger 
Coulomb strength (smaller energy gap) leads to stronger 
electron-electron correlation and consequently the elec- 
trons are more clearly separated from each other. It re- 
sults into a high probability for finding the electrons to 
occupy higher energy states. That also means that the 
system transfers to a configuration with larger Sz and Lz 
at smaller applied field. 

We also found that the ground-state energy is sensitive 
to the presence of the magnetic field. In zero magnetic 
field, the ground state receives contributions from many 
different configurations having z-projection of the mag- 
netic ion, Mz, from —5/2 to 5/2. When a magnetic field 
is applied, the ground state favors states with projection 
of the spin of the magnetic ion down and the states with 
Mz ~ —5/2 give the main contribution to the ground 
state. 

Fig. IHb) shows the average of the three quantities 
Mz , Sz , and Lz of the three-electron quantum dot as a 
function of the magnetic field for a Coulomb interaction 
strength Ac = 0.5. We realize that with increasing Ac 
the < Lz > {< Sz >) exhibits jumps at smaller critical 
flc (compare green dash-dot-dot curve in Fig.[T]). 
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FIG. 7: (Color online) Exchange interaction of a five-electron 
quantum dot with the magnetic ion located at several posi- 
tions and the Coulomb strength Ac = 0.2. The horizontal 
line separates the plane into ferromagnetic (FM) and anti- 
ferromagnetic (AFM) regions. The inset shows Sz alone vs. 
magnetic field. 
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FIG. 8; (Color online) The magnetic ion position dependence 
of the exchange interaction of a three-electron quantum dot 
for Coulomb strengths Ac — 0.2, 0.5, and 2 at fic = 0.4. The 
horizontal dotted line at A^^'^ = separates the plane into 
ferromagnetic (FM) and antiferromagnetic (AFM) regions. 



B. Antiferromagnetic coupling 

Now we direct our attention to the very small mag- 
netic field behavior. There exists a very small region of 
the magnetic field where the total spin of electrons and 
the total spin of the magnetic ion are oriented paral- 
lel, we found this earlier in Ref.l- for the zero magnetic 



FIG. 9; (Color online) Phase diagram for the ferromagnetic- 
antiferromagnetic transition of a three-electron quantum dot 
with the magnetic ion located at (0, 0) (black squares referred 
to (1)) and (2Zo = 5.29A, 0) (red circles referred to (2)). 



field case. These results are now extended to nonzero 
magnetic field. This is made more clear in the upper in- 
set of Fig. [3] where the crossing point of the two terms: 
the Zeeman energy of the magnetic ion (A^'^") and the 
exchange interaction (A|^'^), occurs at fic — 0.01 (con- 
verted Xo B O.IT for the considered system). This 
ferromagnetic couphng extends further, up to = 0.04 
(see the lower inset of Fig. [3]). 

It is worth noticing that this ferromagnetic coupling is 
extended to a much larger magnetic field range (up to 
f2c = 2.3) if we move the magnetic ion to the center of 
the quantum dot (see Fig. |4]). This can be understood as 
follows. When the magnetic ion is located at the center 
of the dot and the magnetic field is very small the ab- 
solute value of A^^" always dominates over A™"^. This 
is opposite to the case when the magnetic ion is located 
at (0.5/o,0). Recall that in Rcf.^^ we found for zero field 
that the exchange Zeeman energy is minimum when the 
magnetic ion is at the center of the quantum dot and ap- 
proximately zero at positions very far from the center of 
the quantum dot. 

Fig. [4] tells us that the magnetic field where the an- 
tiferromagnetic coupling between electron and the mag- 
netic ion starts depends on the position of the magnetic 
ion in the quantum dot. The system with the magnetic 
impurity located at (ZqjO) exhibits an antiferromagnetic 
coupling for fie > 0.1, that is larger than fie = 0.05 in 
case of (0.5/o, 0). 

We have discussed the appearance of antiferromag- 
netism in a three-electron quantum dot. Now we go back 
to the two simpler cases with the number of electrons 
A'e = 1,2 (see Fig. [S]). Let us first discuss the results for 
A'e = 1 as given in Fig.[5Ja). The antiferromagnetic cou- 
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FIG. 10: (Color online) Phase diagram for ferro- 
antiferromagnetic transition of a three-electron quantum dot 
in an applied field of Qc = 0.4 (black squares referred to (1)) 
and Sic = 1.9 (red circles referred to (2)). 



pling between the electron and the magnetic ion starts 
at smaller magnetic field as the magnetic ion is moved. 
This is different from the previous results for A^e = 3. 
The reason is as follows: for the quantum dot with a 
single electron, the electron tends to accommodate per- 
manently the s-shell with < >= (see the inset of 
Fig. [Sl^a)) in the ground state while the exchange param- 
eter in the s shell (Jgs) is found to be maximum right at 
the center of the quantum dot. Moving the magnetic ion 
away from the center of the dot, this Jss is found to be 
smaller and as a consequence the exchange electron- Afn 
interaction becomes smaller than the electron Zeeman en- 
ergy, leading to an antiferromagnetic coupling at smaller 
magnetic field. 

The story for = 2 electrons (see Fig.[5l^b)) is now in- 
teresting since the two electrons accommodate the s shell 
with spins antiparallel making the total spin of the elec- 
tron zero in the ground state with almost unit probabil- 
ity. This leads to zero contribution to the first term writ- 
ten in the last line in [2] for diagonal elements. Therefore, 
the main contribution (even very small) to the exchange 
energy is now expected to come from coupling with con- 
figurations where one of the electrons (spin down) stays 
in the s level and the other occupies higher level (see 
the schematic diagram in Fig. EJ^b)). In this diagram, 
the magnetic ion is assumed to be located at the center 
of the dot with spin down (—5/2). The coupling of the 
electron (spin up) in the s orbital with an electron from 
either of the p shell is zero. The only non-zero coupling is 
with an electron with the quantum numbers (1,0) of the 
d shell (as shown in the diagram and this quantum state 
would change if the ion is located away from the center of 
the dot) with the amount of about — 10~^. This picture 




FIG. 11: (Color online) Radial density scaled by Iq calculated 
for a three-electron quantum dot with the magnetic ion (white 
cross) located at (0.5/o,0) in (a), (b) and (0.8Zo,0) in (c) when 
the Coulomb strength Ac = 0.5 in a magnetic field Qc ~ 6.4 
(see (a)) and 0,^ — 7 (see (b) and (c)). (d) is the magnetic 
ion position dependence of — Af"'^ with the magenta circles 
indicating the respective calculated positions of the magnetic 
ion in (a), (b), and (c). Red color region is for highest value 
of radial density. 



remains valid until the magnetic field is high enough to 
excite one electron from the s shell to a higher quantum 
state forming the ground state with two up spins anti- 
ferromagnetically coupling with the magnetic ion. For 
smaller Coulomb interaction strength, the antiferromag- 
netic behavior occurs at larger magnetic field since the 
two electrons repel each other less and consequently they 
stay longer antiparallel in the s shell. 

From the N^. = 2 result we may ask a question: 
whether a four-electron quantum dot has similar proper- 
ties since the number of electrons in both cases are even 
and one may have a situation where the total spin of the 
electrons is zero. Indeed, if the magnetic ion is located 
at the center of the quantum dot even though for = 4 
the outer shell is half-filled, this is possible as illustrated 
in Fig. [S] where the antiferromagnetic coupling occurs at 
ilc = 3.87 at which the total spin of the electrons reaches 
the maximum value < 5*2 >= 2. In this case, the first 
two electrons will occupy the s shell and the remaining 
two will occupy two of the five orbitals of the p and d 
shells. This picture holds at small magnetic field. How- 
ever, there is a big difference in the exchange energy as 
compared to the previous case of iVe = 2 when the ion is 
shifted away from the center of the dot, e.g. in this plot 
at (0.5Zo, 0). The exchange energy is much larger than the 
result obtained for Rmu = (0, 0) because when the ion is 
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out of the center of the quantum dot, the two remaining 
electrons at higher orbitals have a non-zero contribution 
in the diagonal exchange elements dominating the ex- 
change energy of the ground state. This is the reason 
why the antiferromagnetic transition occurs at smaller 
magnetic field {ilc = 2.5) as compared to the case when 
the ion is at the center of the dot (17c = 3.87) although 
the pictures of the and Lz transition in these cases 
are similar. 

To complete the picture for few-electron quantum dot 
system, we will discuss the AFM behavior for the system 
with the highest number of electrons, Ne = 5, where we 
were able to obtain accurate numerical results. We focus 
on the small magnetic field region (see Fig. [7|). For the 
magnetic ion in the center of the dot, the FM coupling 
is dominant in the shown magnetic field region because 
the diagonal exchange matrix elements dominate over the 
Zeeman energies of the electrons and of the magnetic 
ion. This is different for the cases with the magnetic 
ion displaced a bit from the center of the quantum dot. 
The FM-AFM occurs at Qc = 0.07 and fie = 0.12 for 
Rmu = (0.5^0,0) and Rmu = (^o,0), respectively. It is 
similar to the cases for the system with = 1,3 due 
to the zero coupling between the orbitals from the p, d 
shell with the s orbital. To observe the AFM behavior 
for the system with the magnetic ion located at the cen- 
ter of the dot where the diagonal exchange elements are 
almost zero, it is crucial to include enough quantum or- 
bitals (that rapidly increases the size of the Hamiltonian 
matrix resulting in very time consuming calculations) so 
that one allows the electrons to jump to higher energy 
levels and having parallel spins as previously shown for 
the case = A (see Fig. [B]). In that case, the four- 
electron system exhibits an anti-ferromagnetic coupling 
with the magnetic ion at the magnetic field where the 
total z-projection of the spin is maximum Sz — 2. The 
system is strongly polarized. For the case = 5, up to 
ric = 0.2, the total ^2 = 0.5 and the total Lz = I. The 
inset in Fig. [7] supports the AFM behavior for the out of 
center AIn as obtained in the main plot. 



C. Phase diagram for the 
ferromagnetic-antiferromagnetic transition 

Now we change the Coulomb interaction strength and 
explore the magnetic behavior as function of the position 
of the magnetic ion. From Fig. [51 it is clear that when 
reducing the Coulomb interaction the system undergoes 
a ferromagnetic to antiferromagnetic transition at grad- 
ually larger magnetic fields for Mn^+ positions that are 
closer to the center of the dot. We see that — A^^'^ has 
a peak structure with a maximum at some specific posi- 
tion of the magnetic ion, e.g. see the peak for the case 
Ac = 0.2 (the blue full circles). However, it is certain 
that at high magnetic field, the system is always antifer- 
romagnetic. 

The FM-AFM phase diagram for a three-electron 



quantum dot in (Ac^^c) space is shown in Fig. [9] for 
two different positions of the magnetic ion. When the 
magnetic ion is in the center of the quantum dot (black 
curve with squares in Fig. [5]) the critical magnetic field 
increases as the Coulomb interaction strength decreases. 
The reason is that increasing the Coulomb interaction 
helps the electrons to approach closer the magnetic ion 
and therefore the critical magnetic field for the system to 
transit to the antiferromagnetic phase decreases. 

Now we move the magnetic ion away from the center 
of the quantum dot and we obtain the phase diagram 
as shown by the red curve (with solid circles) in Fig. (S] 
For Ac < 0.4, the stability of the FM phase with re- 
spect to an applied magnetic field is strongly reduced 
and a small magnetic field turns the three-electron sys- 
tem into the AFM phase. Notice that for sufficient strong 
electron-electron interaction (i.e. Ac > 0.4) we obtain 
practically the same FM-AFM phase diagram as for the 
case the Mn-ion is located in the center of the quantum 
dot. A remarkable re-entrant behavior is found in the 
region 0.3 < Ac < 2 and 0.9 < ilc < 2.5 where with 
increasing Ac we go from an antiferromagnetic to ferro- 
magnetic and back to antiferromagnetic phase. This un- 
usual behavior is understood as follows. As the impurity 
is moved away from the center of the quantum dot the ex- 
change matrix will have many nonzero off-diagonal terms 
that leads to a smaller FM-AFM critical transition mag- 
netic field. Now let us turn our attention to the region 
Ac < 0.4. For very small Coulomb interaction strength 
the electrons will repel each other only weakly and are 
therefore pulled towards the magnetic ion (the nonzero 
exchange matrix elements increase strongly) resulting in 
a very small FM-AFM magnetic field. For Ac > 0.4, the 
electrons become more strongly correlated and the criti- 
cal field stays about Ailc ~ 0-02 0.07 from the result 
for Rmu = 0. If one moves the ion further and further 
away from the center, the Ac ~ 0.4 transition line moves 
to larger Ac values. For example, for Ac = 0.5 and the 
magnetic ion located at (0.5^o,0) the FM-AFM critical 
transition occurs at flc — 0.08, which is much smaller 
than 2.21 found for Rmu = (0.2^0,0). 

The dependence of the ferromagnetic- 
antiferromagnetic transition of a three-electron quantum 
dot system on the position of the magnetic ion is 
summarized in the phase diagram shown in Fig. [TO] for 
two different magnetic fields fie = 0.4,1.9. We can 
predict that with slightly larger (smaller) magnetic field 
the slope of the curve will be larger (smaller). From 
Fig. [51 we already learned that the FM-AFM transition 
magnetic field is largest for the ion in the center of the 
quantum dot as is also seen in Fig. [TOl The re-entrant 
behavior of the AFM phase as function of Ac is found 
for small Rmu/Io values, i.e. when the Mn-ion is not 
too far from the center of the quantum dot, in case 
the magnetic field is not too small. The critical point 
(^Mn/^Oj Ac) = (0. 457, 0.656) for = 1-9 moves down 
(up) with increasing (decreasing) magnetic field. 
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FIG. 12: (Color online) Spin up-up pair correlation functions 
scaled by Iq calculated for the quantum dot of Fig. 1111 where 
one electron (orange arrow) is pinned at one of the three most- 
probable positions of the electrons as obtained in Fig. 1111 It 
is at (-0.44^0,0) in (a) for Fig.IlUa), (-0.63^0,0) in (b) for 
Fig. [nib), and (-0.52/o,0) in (c) for Fig.fTlTc'). respectively, 
(b') The same correlation function for the case the position 
of the fixed electron is closer to the ion (at (— 0.3Zo,0)) as 
compared to (b). The position of the magnetic ion is indicated 
by the white cross. 



D. Density and correlation 

In high magnetic field, the magnetic ion tends to at- 
tract electrons because they are oppositely polarized. Be- 
cause the exchange interaction is small as compared to 
the Coulomb interaction, the electrons and magnetic ion 
are arranged in such a way that the electrons repel each 
other and also try to be as close to the magnetic ion as 
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FIG. 13: (Color online) Addition energy for different number 
of electrons. The magnetic ion is located at (0.41Zo, 0) and the 
Coulomb interaction strength is Ac ~ 0.2. 

possible. This picture holds above the FM-AFM critical 
magnetic field. 

To show this behavior explicitly, we studied the radial 
density and the radial pair correlation functions. Their 
respective operators are defined as: 

i=l 

and 

(11) 

We plot in Fig. [11] the radial density of a three electron 
quantum dot that is polarized in high magnetic field for 
the case that the Coulomb strength is Ac — 0.5 and 
the magnetic ion is located at two different positions 
for two magnetic fields. The electrons and the magnetic 
ion are antiferromagnetically coupled. The strength 
of that coupling can be seen from Fig. [TlTd) in which 
we plot the magnetic ion's position dependence of the 
exchange energy at two magnetic fields i7c = 6.4, and 
7. Those magnetic fields are typical in the sense that 
the exchange term is found to be very large (Jlc = 6.4) 
or the correlation between the electrons very high 
(fic = 7). Density plots are shown for Run at (0.5Zo,0) 
and (0.8^0,0). We observe three distinct peaks of max- 
imum probability. They are found at: (— 0.44Zo,0), 
(0.22,0.44)^0, and (0.22,-0.44)^0 in Fig. [TTIa); 
(-0.63^0,0), (0.26,0.63)^0, and (0.26,-0.63)^0 in 
Fig-Elb); (-0.52^0, 0), (0.26, 0.52)/o, and (0.26, -0.52)^0 
in Fig. [TUc). These figures show clearly the interplay 
effect where the three electrons on the one hand try 
to be close to the magnetic ion and on the other hand 
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repel each other via the Coulomb potential energy. It 
results in the merging of the radial density such that 
the higher the exchange energy the larger the merging 
of the local maxima in the electron density and the 
smaller the correlations. Fig. [TlT d) gives an idea about 
the variation of — Af""^ with the position of the magnetic 
ion and reaches a maximum at (O.SZq, 0) for fie = 7. 
In Fig. [TTT c). the three electrons are less attracted to 
the magnetic ion via the antiferromagnetic coupling as 
compared to that in Fig. [TIT b). This is due to the fact 
that the — A^'^ for the case shown in Fig. fllTb) is larger 
than that in Fig. [TTT c). The electrons are therefore 
found more correlated in the latter case presented by the 
extended red region in Fig. [TlT c). Thereby, correlation 
between electrons in Fig. [TlT c) is expected to be the 
highest and in Fig. [TTT a) the smallest. 

The position of the magnetic ion affects the ground- 
state property as is made clear in Fig. [T^l We fix the spin 
state and the position of one electron (indicated by the 
orange arrow) and the position of the magnetic ion (the 
white cross). The magnetic field is such that flc — 6.4 in 
Fig- flZT a) and flc = 7 in the others. It also reflects the 
fact that the system in Fig. [TIT a) exhibits the smallest 
correlation as compared to the other two. This illustrates 
the point raised above about the density. At the magnetic 
field = 7, the electrons are strongly polarized result- 
ing in the red regions of the up-up spin pair correlation 
function that tends to surround the magnetic ion. We see 
that the three electrons arc most likely to localize around 
some specific positions defining a triangle with the three 
electrons at the three vertices, while they are attracted 
to the magnetic ion. When we locate one electron at a 
position closer to the magnetic ion, see Fig. [TWb'). the 
two peaks decrease in amplitude as compared to those in 
Fig.mb). 



E. Addition energy 

The addition energy (often called the chemical po- 
tential) is defined as the increase of the energy of the 
quantum dot system when an electron is added: iin^ — 
Easi^e) — EGsi^e — 1)- This quantity can be measured 
experimentally and is plotted in Fig. [T2] as function of 
the magnetic field. 

There are several cusps appearing in the addition en- 
ergy curves as a consequence of changes in the ground 
state. These changes are due to variations in the z- 
projection of the total spin of the electrons and/or the 
z-projection of the total angular momentum of the sys- 
tem when the magnetic field increases beyond some spe- 
cific values. The presence of the magnetic ion leads to 
more cusps and the position of these cusps is also influ- 
enced by the number of electrons and the position of the 
magnetic ion. The cusps are from either of the two sys- 
tems in the study. For instance, the green triangles in 
Fig. [HI are for ^3 = EasiNe = 3) - EasiNe = 2) has 
two cusps at ilc — 2.6 and ilc = 3.4. The cusp at the 




-2-10 1 2 

FIG. 14: (Color online) Contour plot for the phase of the re- 
duced wave function of a three-electron quantum dot in mag- 
netic field f2c = 11 with the magnetic ion (the white cross) lo- 
cated at (0.5Zo, 0) in (a),(b) and at (O.SZo, 0) in (c) for Coulomb 
strength Ac ~ 0.5. Two fixed electrons (indicated by the two 
magenta circles) are located at the two peaks appearing in 
the radial density: (—0.63^0,0) and (0.26, 0.63)/o in (a), (c); 
and (-0. 63^0,0) and (0.26, -0.63)^0 in (b). 



point flc = 2.6 comes from the change in the configu- 
ration of the average of the total z-projection spin and 
the total z-projection of angular momentum (Sz,Lz) of 
the two-electrons in the quantum dot from (0,0) to (1,1). 
While the other cusp fie = 3.4 comes from the change of 
the phase of the three-electron quantum dot from (0.5, 1) 
to (1.5,3). It is similar to the case for /i4 (the blue left 
pointing triangles) where the cusp appears at Qc = 3.4. 
At this point, we observed a change from configuration 
(0,2) to (1,3). The remaining one, fie = 3.9, is from the 
four-electron case when its configuration changes from 
(1,3) to (2,6). 



F. Vortex structure: many-body correlations 

Another way to obtain information on the correlations 
that are present in the many-particle wave function is to 
investigate the vortex structure. At a vortex the many- 
body wave function is zero and is characterized by a 
change of phase of 2t: when we go around this point. 

The zeros of the wave function are similar to flux 
quanta when e.g. the wave function corresponds to the 
order parameter in a superconductor. The fixed electrons 
and the zero of the wave function follows closely the dis- 
placed electron and one may say that the electron plus its 
zero form a composite-fermion object. The composite- 
fermioni^iiiii^ (and references therein) is a collective 
quasi-particle that consists of one electron bound to an 
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FIG. 15: (Color online) The energy spectrum of the first 30 
levels of the three (the top), four (the middle), and five (the 
bottom) electron quantum dot with the magnetic impurity 
located at (0.5io,0). The Coulomb strength is Ac = 0.5 with 
magnetic fields flc = 0.0 (the blue squares) and Qc = 0.4 (the 
magenta circles) . The inset in the top figure shows the average 
of Mz (the dark-blue triangles) and Sz (the red triangles) as 
a function of Q,c of the sixth energy level. 
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FIG. 16: (Color online) The magnetic field dependence of the 
energy spectrum (the first 120 levels are shown) of the one- 
electron quantum dot with the magnetic impurity located at 
(0.5/o,0) and Ac = 0.5 (the inset is a zoom of the first 24 
levels at low magnetic fields and the inset in the inset shows 
the first 12 levels, it is a zoom of the rectangular region. 



even number of vortices (flux quanta). The composite- 
fermion concept introduces a new type of quasi-particle 
that is used to understand the fractional quantum Hall 
effect in terms of the integer quantum Hall effect of these 
composite fermions. 

To obtain the zeros of the wave function of the sys- 
tem with A^e electrons, we fix A^e ~ 1 electrons at some 
positions inside the quantum dot and leave the remain- 
ing one free. The resulting reduced wave function gives 
the probability to find the remaining electrons at differ- 
ent positions in the quantum dot and the zeros' of this 
function are those points where the phase of the wave 
function changes by 27r. As an example, we investigate 
the situation of a three-electron quantum dot. 

Figs. [13] shows the vortex pictures of a three-electron 
quantum dot containing a magnetic ion located at posi- 
tions (see the white cross) that are identical to its posi- 
tions in Fig. lllf a). Two, among three, electrons are fixed 
at the respective peaks in the electron density. Red and 
black regions are referring to the highest (27r) and lowest 
phase 0, respectively. Those plots show that there are 
always two vortices near the pinned electrons' positions. 
For example, the number of vortices pinned to each elec- 
tron in the case iVe = 3 at fie = 11 is 2 describing the 
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FIG. 17: (Color online) The same as Fig.[T6]but now for the 
two-electron quantum dot. The inset in the inset is a zoom 
of levels 7th to 24th for magnetic field close to zero. 

system at filling factor v = ^2<l~>^ ~ ^Z"^- Notice also 
that one of the vortices appears to be pinned at a position 
very close to the Mn ion. 

We realize that moving the magnetic impurity to a dif- 
ferent position changes the relative positions of the vor- 
tices that are pinned to the electrons with respect to one 
another, as shown in Figs. [T^a') and (c). As the elec- 
trons are antiferromagnetically coupling to the magnetic 
ion this kind of movement consequently depends on the 
position of the magnetic ion. 

In the case Ac = 0.5, we found that the average of the 
maximum z-projection of the total angular momentum is 
< Lz >= 9 and the two vortices appear at the external 
field fie = 3.0. Apparently, the larger Ac the smaller 
ric for which the first two vortices appear at the pinned 
electrons. 



IV. ENERGY SPECTRUM 

In the presence of an external magnetic field, the many- 
fold degeneracy of the energy spectrum of the system is 
lifted. Fig. [T5f a) illustrates that point for the case of 
three electrons. In the absence of the interaction between 
the electrons and the magnetic ion and in the absence of 
a magnetic field (blue squares), the energy spectrum is 
7-fold degenerate for the first 7 lowest energy levels, the 



next level is then 5-fold degenerate, and the next 7-fold 
degenerate, and so on. The origin of this was explained in 
Ref.— and is due to the coupling of the electrons and the 
magnetic ion. When the magnetic field is different from 
zero, see red circles in Fig. [TSl the degeneracy is lifted. 
In the inset of Fig. [15] we plot < > (the magenta 
triangles) and < Sz > (the dark-blue ones) as a function 
of magnetic field for the sixth level. The average of < 
Mz > and < Sz > change abruptly as compared to those 
found for the ground-state energy, e.g. see Fig. [2l Mz 
and < 5*2 > of the sixth state jump between two different 
values, e.g. —1.5 and —2.5 for < Mz > and 0.5 and 1.5 
for < 5*2 > as function of the field. This is a consequence 
of anti-crossings of energy levels as will be apparent later. 
The result for four- and five-electron quantum dots are 
also shown in Fig. [151 We see the degeneracy of 8, 6, 4, 
and 12 for the first 30 levels in the case = 4 and of 7, 
5, 7, 5, 5, and level 30 has the same degeneracy with the 
next energy level beyond the first 30 for the case = 5 
in B = Tesla. 

To have a clearer picture of the energy spectrum of 
the quantum dot system we plot in Figs. [THl [T71 and [TH| 
the magnetic field dependence of the first 120 energy lev- 
els for Ne = 1,2,3, respectively. The spectra at small 
magnetic fields is enlarged (see insets) to see the Zeeman 
splitting and the nearly-linear behavior of the energy lev- 
els. Remember that this is due to the coupling of the 
electron spins with the magnetic ion spin. For iVg = 1 
the first two levels for B = are 7- and 5-fold degenerate 
(7-fold degenerate is due to the ferromagnetic coupling of 
the s-shell electron spin 1/2 with the magnetic ion spin 
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5/2 and 5-fold degenerate of that electron now with spin 
— 1/2 to the magnetic ion with spin 5/2), respectively. A 
closer inspection (see Fig. [16]) tells us that these 12-levels 
are exchange split into two bundles of 7- and 5-levels (the 
inset in the inset in Fig. [T6)) . Notice that there is a first 
large energy gap at very small fields between the first 12 
levels and the next 24 ones as seen in Fig. [11] while that 
kind of gap appears between the first 6 and the next 36 for 
TVe = 2 (Fig.[T7]). For TVe = 3 (Fig. [18]), this kind of gap 
appears after the first 24 energy levels. The origin is the 
coupling of the third electron, which can reside at either 
two states of the p shell while the s shell is already fully 
filled, with the magnetic ion with 6 z-components of the 
spin at very small fields, i.e. the intra-shell (p) exchange 
interaction. For Ne ~ 2 the electron ground state corre- 
sponds with a filled s-shell, i.e. < Sz >— 0, and therefor 
for S = only a 6-fold degeneracy, as shown in Fig. [TTl 
is found due to the z-component of the Mrt-spin. The 
next level is 8-fold degenerate at B = (8 come from the 
ferromagnetic coupling of the two-electron system with 
total spin 1 to the magnetic ion with spin 5/2) (see the 
inset in the inset in Fig. [T7)l . etc. 

With increasing magnetic field, we see that for — 1 
there is periodically an opening of energy gaps in the 
spectrum. Similar energy gaps have been found earlier 
(as an example see e.g Ref.^^) for a quantum dot with- 
out a magnetic impurity and are a consequence of the 
electron with two-fold spin degeneracy filling the equally- 
gaped-energy single-particle quantum states with differ- 
ent sets of the radial and angular quantum numbers. No- 
tice that for Ne = 2, 3, these gaps have disappeared. 

The spectra exhibit a lot of crossings and anti- 
crossings, the number of them has increased as compared 
to the quantum dot case without a magnetic ion because 
of the Zeeman splitting of the Mn-spin. When the ap- 
plied field increases the gaps in the spectrum of A^e = 1 
are still open and appear more often than in the cases of 
Ne = 2,3. Once again, we see a lot of cusps in the energy 
levels and that reminds us to abrupt changes in the con- 
figuration of the system with magnetic field as discussed 
before for the ground state. 



THERMODYNAMIC PROPERTIES 



A. Magnetization and susceptibility 
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FIG. 19: (Color online) Magnetization of the ground state 
of the two- and three-electron quantum dot with the mag- 
netic ion located at (0.5/o, 0) for three values of the Coulomb 
strength Ac ~ 0.2, 0.5, and 1.1. 



also has a peak at flc = 3.3. The same thing happens at 
ilc = 1.4,4.1, and 6.8 for Ac = 1.1 in the magnetization 
and the susceptibility. 

For non-zero temperature, the temperature depen- 
dence of the magnetization and susceptibility is defined 
by: M{T) = -d{E{T))/dB, x{T) = dM{T)/dB, respec- 
tively. The statistical average < E{T) > is calculated as: 



We first calculate the magnetization and susceptibility 
of the system: M = -dEcs/dB and x = dM/dB at 
zero temperature. 

The magnetization of a quantum dot with the magnetic 
ion located at (0.5Zo, 0) with A'e = 2, 3 electrons is plotted 
in Fig. [TH] We see several jumps that are a consequence 
of changes in the ground state, e.g. changes in < > 
(see previous section) . For example, the magnetization of 
the three-electron quantum dot as plotted in Fig. [TWb) 
for the case Ac = 0.2 and the magnetic ion at (0.5Zo,0) 
has a step at flc = 3.3. Consequently, the susceptibility 



{EiX,T, Rmu)) 



(12) 



where the sum is over the energy levels as displayed in 
e.g. Fig.fTSl 

These quantities are explored in Figs. [20] for A^e = 3 
and a few different temperatures (including the zero tem- 
perature case). With increasing temperature the jumps 
become smoother. A very low magnetic field peak shows 
up because for T 7^ we have Af w at fij. = 0. 
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FIG. 20: (Color online) Temperature dependence of the mag- 
netization and susceptibility for a three-electron quantum dot 
with the magnetic ion located at (0.5/o,0) and Ac = 0.2. 



B. Heat capacity 

An important quantity that is related to the storage of 
energy is the heat capacity: 



Cv (A, T, Rmu) 



d{E{X,T, Rmu)) 
dT 



(13) 



The heat capacity is investigated as a function of the 
Coulomb strength A, temperature T, magnetic field, and 
the position of the impurity Rmu- 

We plot in Fig. [^1] the specific heat for two values of 
the magnetic field, i.e. f2c = 0.5 and f2c — 7.5, and 
three typical positions of the magnetic ion. For weak 
fields, the three electrons start to polarize and we see 
that the position of the main peak moves towards higher 
temperature as the magnetic ion is moved away but not 
too far from the center of the dot. For the high magnetic 
field case the three electrons are strongly polarized and 
we see a different behavior in the shift of the main peak. 
This results from the change of the statistical average of 
the energy levels at different fields. 



-I— • 

o 

CD 
Q. 
CD 
O 



" 1 



CD 
0) 
X 







1 1 


- 


RMn(0-5lo.O) 






RMndo-O)^ 






Ly Q =0.5 > 

— ^ c\ 








Rm (0.0) 






RMn(0-5io>o) : 


I 





1 . 1 



20 



40 60 
T(K) 



80 



100 



FIG. 21: (Color online) The heat capacity vs temperature of 
a three-electron quantum dot with the magnetic ion located 
at three different positions, the Coulomb interaction strength 
is 1.0 and the magnetic field Q.c ~ 0.5 (the three lower curves) 
and 7.5 (the three upper curves). 
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FIG. 22: (Color online) The heat capacity vs magnetic field 
of a three-electron quantum dot with the magnetic ion located 
at (0.41Zo,0), Ac ~ 0.2, at several temperatures: 23, 116, 174 
and 209 K. 



Now we examine the behavior of the heat capac- 
ity at a specific temperature as function of magnetic 
field. Figs. [22l [23l and [24] are the plots of the magnetic 
field dependence of the heat capacity of three (at two 
Coulomb interaction regimes) and four-electron quan- 
tum dots at some specific temperatures and two different 
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FIG. 23: (Color online) The heat capacity vs magnetic field 
of a three-electron quantum dot with the magnetic ion located 
at (0.5Zo,0), Ac ~ 0.5, and the temperatures are 5.8K and 
U.6K. 
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FIG. 24: (Color online) The heat capacity vs magnetic field 
of a four-electron quantum dot with the magnetic ion located 
at (0.41/o,0), Ac = 0.2 plotted for several temperatures. 



VI. DISCUSSIONS. 



Ac = 0.2,0.5. The peak at small magnetic fields broad- 
ens and moves to higher fields with increasing tempera- 
ture. The heat capacity exhibits a number of peaks and a 
clear minimum around e.g. Wc = 3.4 as shown in Figs. 1221 
Remember that this field corresponds to a cusp in the 
energy versus magnetic field behavior as discussed previ- 
ously in subsection III. A. At very low temperatures, this 
cusp still affects the heat capacity through the sharpness 
of the minimum as shown in the figure and this gradually 
becomes small at high temperatures. In Fig. [211 we see a 
very interesting behavior of the heat capacity at Qc ~ 2.3: 
the single peak becomes a double peak as temperature in- 
creases from T = 5.8 to T = 11.6K. This is due to the 
cusps now occurring around this field in the low-energy 
levels of the spectrum of the three-electron quantum dot 
system as observed in Fig. [THl Besides, the structure of 
the heat capacity is more complex (more peaks) with in- 
creasing Ac. This is made clear if one looks back to the 
previous discussion related to Figs. [1] and [U^a). 

For the case Ne = 4, the heat capacity exhibits more 
peaks as compared to the case = 3 and the behav- 
ior of the peaks with increasing temperature is also very 
different. Temperature affects the heat capacity of the 
system in the sense that it increases the peak values and 
separates them in magnetic field. 

The Coulomb interaction strength changes the struc- 
ture of the magnetic field dependence of the heat capac- 
ity and is illusive through Figs. [22] and [231 The peak 
of the heat capacity for the case with smaller Coulomb 
interaction strength appears at higher magnetic field as 
compared to the case with larger one. 



Due to the presence of the magnetic ion (and electron- 
electron interaction) , electrons in the ground state do not 
always completely polarize in the presence of an external 
magnetic field. The configurations are mixed consisting 
of electrons having spins up and down. But for very 
large magnetic field, the magnetic ion tends to pull the 
electrons closer to the ion forming a ring-like electron 
density profile. These are the consequences of the inter- 
play of several effects such as the Zeeman effect (on the 
electrons' and the magnetic ion's spins), the Coulomb 
repulsion, the spin exchange interaction, etc. This com- 
petition results in a crossover from ferromagnetic- to an- 
tiferromagnetic coupling between the electrons and the 
magnetic ion at some specific magnetic field. Interest- 
ingly, this effect is observed to appear at higher magnetic 
field when we move the magnetic ion further from the 
origin of the quantum dot. A re-entrant behavior of the 
FM-AFM transition is found as function of the Coulomb 
interaction strength when the magnetic ion is moved out 
(but not too far) from the center of the quantum dot. 

The energy levels exhibit cusps as function of the 
magnetic field which correspond to changes of the con- 
figuration of the system as expressed by the values of 
(< Sz >,< Lz >). These cusps move to lower mag- 
netic field with increasing Coulomb interaction strength. 
The number of cusps increases with increasing number of 
electrons. These cusps show up in the addition energy. 

The transformation of the electron system to those of 
composite fermions is studied. In high magnetic fields, 
the electrons attach an even number of quantized vor- 
tices which we made clear by examining the many-body 
ground-state wave function in the presence of a magnetic 
ion. Unlike the case without a magnetic ion where all 
the vortices are tightly bound to the electrons, when we 
fix the electrons at different positions the system of vor- 
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tices stays pinned to the electrons and moves with the 
electrons but the relative positions of the vortices are 
modified. 

The contribution of the local Zccman splitting energy 
to the total energy of the system in large external fields 
is very small as compared to the contributions from the 
other parts. However, a slight movement of the position 
of the magnetic ion inside the quantum dot affects the 
result, shghtly. 

With increasing applied magnetic field, each time the 
system jumps to a different (< >,< Sz >) config- 
uration leads to the appearance of a peak in the ther- 
modynamic quantities as e.g. the susceptibility and the 



heat capacity. In the presence of the magnetic ion, the 
structure of peaks in the heat capacity changes with the 
position of the magnetic ion. As temperature increases, 
these peaks split into two peaks and become smoother. 



VII. ACKNOWLEDGMENTS 

This work was supported by FWO-Vl (Flemish Science 
Foundation), the EU Network of Excellence: SANDiE, 
and the Belgian Science Policy (lAP). 



* Electronic address: lnga.nguyen@ua.ac.bel 

^ Electronic address: francois.peeters@ua.ac. be] 
^ J. K. Furdyna, J. Appl. Phys. 64, R29 (1988). 
^ D. L. Klein, R. Roth, A. K. L. Lim, A. P. Alivisatos, and 

P. L. McEuen, Nature (London) 389, 699 (1997). 
^ Y. F. Chen, J. H. Huang, W. N. Lee, T. S. Chin, R. T. 

Huang, F. R. Chen, J. J. Kai, and H. C. Ku, Appl. Phys. 

Lett. 90, 022505 (2007). 
^ K. Chang, J. B. Xia, and F. M. Peeters, Appl. Phys. Lett. 

82, 2661 (2003). 
^ A. O. Govorov, Phys. Rev. B 72, 075358 (2005). 
® Paul I. Archer, Steven A. Santangelo, and Daniel R. 

Gamelin, Nano Lett. 7, 1037 (2007). 

J. Fernandez-Rossier and L. Brey, Phys. Rev. Lett. 93, 
117201 (2004); L. Besombes, Y. Leger, L. Maingauh, D. 
Ferrand, H. Mariette, and J. Cibert, Phys. Rev. Lett. 93, 
207403 (2004); Y. Leger, L. Besombes, L. Maingauh, D. 
Ferrand, and H. Mariette, Phys. Rev. Lett. 95, 047403 
(2005); ibid. Phys. Rev. B 72, 241309(R) (2005); L. Main- 
gault, L. Besombes, Y. Leger, C. Bougerol, and H. Mari- 
ette, Appl. Phys. Lett. 89, 193109 (2006). 

* T. Schmidt, M. Scheibner, L. Worschech, A. Forchel, T. 
Slobodskyy, and L. W. Molenkamp, J. Appl. Phys. 100, 
123109 (2006). 

^ P. Wojnar, J. Suffczyriski, K. Kowalik, A. Golnik, G. Kar- 

czewski, and J. Kossut, Phys. Rev. B 75, 155301 (2007). 
° J. Fernandez-Rossier and Ramon Aguado, Phys. Rev. Lett. 



98, 106805 (2007). 

Y. Leger, L. Besombes, J. Fernandez-Rossier, L. Main- 
gault, and H. Mariette, Phys. Rev. Lett. 97, 107401 (2006); 
L. Maingault, L. Besombes, Y. Leger, H. Mariette, and 
C. Bougerol, Phys. Stat. Sol. (c) 3, 3992 (2006); M. M. 
Glazov, E. L. Ivchenko, L. Besombes, Y. Leger, L. Main- 
gault, and H. Mariette, Phys. Rev. B 75 205313 (2007). 
M. B. Tavernier, E. Anisimovas, F. M. Peeters, B. Szafran, 
J. Adamowski, and S. Bednarek, Phys. Rev. B 68, 205305 
(2003); M. B. Tavernier, E. Anisimovas, and F. M. Peeters, 
Phys. Rev. B 70, 155321 (2004); ibid. Phys. Rev. B 74, 
125305 (2006); T. Stopa, B. Szafran, M. B. Tavernier, and 
F. M. Peeters, Phys. Rev. B 73, 075315 (2006). 
" Shun- Jen Cheng, Phys. Rev. B 72, 235332 (2005). 

Fanyao Qu and Pawel Hawrylak, Phys. Rev. Lett. 95, 
217206 (2005); ibid. 96, 157201 (2006). 
Nga T. T. Nguyen and F. M. Peeters, Phys. Rev. B 76, 
045315 (2007). 

Ramin M. Abolfath, Pawel Hawrylak, and Igor Zutic, 
Phys. Rev. Lett. 98, 207203 (2007). 

J. K. Jain, Phys. Rev. Lett. 63, 199 (1989); ibid. Phys. 
Rev. B 41, 7653 (1990); ibid. 42, 9193(E)(1990). 
^® H. Saarikoski, A. Harju, M. J. Puska, and R. M. Nieminen, 
Phys. Rev. Lett. 93, 116802 (2004). 

S. Tarucha, D. G. Austing, T. Honda, R. J. van der Hage, 
and L. P. Kouwenhoven, Phys. Rev. Lett. 77, 3613 (1996). 



